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FOREWORD 


This  work  constitutes  Part  II  of  a  two  part  final  report  covering  theoretical 
and  experimental  aspects  of  the  general  problem  of  *convective  heat  transfer  with  chem¬ 
ical  reaction.  Part  I  is  directed  toward  the  theoretical  development  of  correlation 

^  • 

formulae  for  the  prediction  of  heat  fluxes  in  high  performance  rocket  motors  and  rela-  * 
ted  systems.  The  research  reported  on  herein  deals  with  experimental  studies  of  heat 
transfer  from  combustion  gases  to  solid  surface^-.  ® 

The  author  has  benefited  greatly  from  frequent  discussion  with  Professor  R. 
Eichhorn  of  Princeton  University  throughout  the  course  of  this  research.  It  is  alro 
a  pleasure  to  thank  Professor  J.®  Fenn  for  his  initial  guidance  and  Ted  Douglass  for 

»  '  S'  ^ 

f)  S' 

his  skill  and  persistence  in  carrying  out  the  subatmospheric  experiments.  ®The  contri.-’ 

*  I®  « 

butions  of  Frank  Kuehner,  Donald  Jost’and^Charlotte  Phillips  in  the  computational  phases 

of  this  research  are  greatfully  acknowledged.  Joseph  Shetz  assisted  with”  the  .assembly  . 

0  # 
of  the  burner  and  mercury  flow  system.  Thanks  are  due  to  Robert  Zeh  for  his  assistance 

with  the  explosion  program  instrumentation.  Finally,  the  author  particularly  wishes  to 

«  * 
thank  Dr.  H.  F.  Calcote  of  this  laboratory  and  Mr.  Everett  Stephens  and  Mr.  Erich 

0 

Soehngen,  “Of  the  Aeronautical  Research  Laboratory,  Office  of  Aerospace  Research, 

'  *  0 

for  their  technical  interest  and  continued  support. 
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NOMENCLATURE 


«  • 

« 


exponent  appearing  in  Eqt(23) 
exponent  appearing  in*Eq.(23)  *  •  • 
specific,heat  at  constant  pressure  ® 

specific  heat  at  constant  volume 

diameter  *  *  • 

•  • 

Pick  diffusivity 

•  •  , 

specific  internal  energy  . 

total  internal  energy  in  vessel  * 

gravitational  acceleration  at  sea  level  % 

Grashof  number,  Eq.{l8)  *  * 

lengtff  over  which  coolant  temperature  rise  is  measured 

mass  flow  rate 

molecular  weight 

summation  index,  Eq.(ll) 

Prandtl  number  for  heat  conduction  =  v/a  ^ 

total  pressure  » 

heat  transfer  rate  per  unit  area 

total  heat  transferred  across  boundary  * 

bomb  radius  * 

universal*gas  constant 

surface  area  of  vessel  walls 


absolute  temperature  ^ 
enclosed  vclume  within  vessel 
thermal  diffusivity 

coefficient  of  thermal  expansion;  1/T  for  a  perfect  gas 
specific  heat  ratio  * 

•  m  ^ 

boundary  layer  thickness  • 

total  hemispheric  emissivity 
thermal  conductivity  of  fluid  or  solid  ' 
kinematic  viscosity  of  fluid  i 

absolute  density 

Stefan-Boltzmann  radiation  constani^ 
nondimensional  time 


Vll  f, 


O  Q  O 


0 


(S'  % 


•  • 

<1>  equivalence  ratio  (fuel  to  oxidizer  ratio  by  volume  or  fiass  normalized  by  the 
corresponding  stoichiometric  ratio)  • 

Subscripts  • 

•  •  •  *  • 

a  apparent  •  ^ 

w 

Hg  pertaining  to  mercury 

f  film  *  ^  » 

k  pertaining  to  species  k 

m  at  the  moment  of  maximum  pressure 

0  at  initial  condition  (usually  vihen  ignition  occurs)  * 

T  partial  derivative  at  constant  temperature 

u  universal 

w  at  the  wall  or  pertaining  to  the  solid 
e  pertaining  to  radiation 

X  pertaining  to  heat  conduction 

p  partial  derivative  at  constant  density 

crit  critical 


Superscripts 

*  normalized  by  the  value  of  the  total  heat  loss  at  t 

t 

Miscellaneous 

n  average  value  of  (  ) 

A  change  in 

5  partial  differential  operator 
d  ordinary  differential  operator 


=  00 


exp(  )  the  exponential  function  e 


(  ) 


Vlll 


*  • 


*  .*» 


t#  • 


e 


f 
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•  *1.  IjJrroDUCTION  •  *.  *  •  .  * 

•  -  •  •  • 

*  *  •  % 

It  has  been  recognized  for  some  time  that  energy  transfer  between*gases 

•  and  solid  surfaces  can  be  affected  by  thermochemical,  processes  ancf,  in  particular, 
the  kinetics  of  these  processes^ This  may  be  viewed  as  either  a  blessing^r 

a  nuisance,  depending  on  one's  point  of  view.  In  either  case,  however,  it  is  nec- 

•  ♦  *  • 

essary  to  understand  #the  magnitude  of  the  effects  to  be  expected,  and  the  way  in 

which  these  effects  will  change  with  changes  in  environmental  conditions.  ,The  de- 

signer  of  those  stationary  power*.stations  relying  ultimately  on  the  energy  content 

•  •  •  • 
of  chemical  fuels  will  always  be  interested  in  finding,  understanding,  and  applying 

methods  of  increasing  the  effitiency  and  rates  of  energy  transfer  from  reactant  or 

•  *  * 

product»gases  tcj|  heat  exchanger  surfaces.  Similarly,  thermochemical  effects'can 

govern  the  choice  of  gaseous  coolants  for  such  devices  a^  nuclear  reactors,  high- 

»  , 

current  electrical  generators,  compact  electromagnets  and  related  equipment.  On 

# 

the  other  ha^d,  in  the*design  of  thermal  protection  systems  for  rocket  motors  or 

•  •  • 

high  speed  aerodynamic  vehicles,  the  problem  of  augmented  heat  transfer  rates  must 

be  anticipated  and  provided  for.  ^  •  * 

♦  Inrfnany  cases,  heat  transfer  considerations  govern  the  design  of  high  per- 

*  i  •  -  • 

formance  propulsion  systems.  Wfth  current  interest  in  more  energetic  chemical  fuels 

* 

•  and  Icfnger  thrust  durations,  the  nozzle  cooling  problem^»^>°  is  correspondingly  in- 

*  • 

tensified.  It  is  for  these  new  propellant  combinations  that  thermochemical  effects 

*  are  likely  to  become  noticeable  since  a  substantial  fraction  of  the  reaction  pro- 

•  •  • 
ducts  at  the  chamber  pressure  and  temperature  are  in  the  form  of  light,  dissociated 

^gases.  Anomalously  large  transfer  rates^have  already  beerf  reported”^  in  such  sys¬ 
tems,  and  thes^,have  been  presumed  to  be  due  *to  tlTe  effects  of  exothermic  chemical 

#  . 

reaction,  both  in  the  gas  phase  and  at  the  gas/solid  interface,  as  well  as  enhanced 

radiation  from  the  chemically  reacting  mixture.  However,  due  to  the  complexity  and 

expense  of  full  scale  rocket  motor  tests,  experimental  results  are  few  in  number 

«  • 
and  very  often  are  of  limited  value  for  purposes  of  generalization.  A  large  price 

has  already  b^en  paid  for  the  experimental  information  currently  available  in  tfie 

« 

literature. 

Taking  the  view  that  a  more  fundamental,  though  perhaps  more  academic, 
set  of  experimental  measurements  would  prove  useful  in  gaining  an  understanding  of 

the  influence  of  chemical  reaction  on  h*eat  transfer,  this  laboratory  embarked  on 

•  *  • 
such  a  program  in  April  of  1959  under  the  support  of  the^ United  States  Air  Force. 

As  a  result,  two  new  experimental  techniques  have  emerged  which  offer  compararf:ive^» 


•  • 
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economy  and  simplicity  and  appear  to  be  capable  of  providing  a  quantitative  under- 

standing  of  the  important  thermochemical  and  radiative  energy  transport  mechanisms. 

»  •  •  o 

In  support  of  these  studiesj  parallel  theoretical  investigations  have  been  carried 

out  on  catalysis  effects  in  flow  systems^*^,  and  convective  heat  transfer  in  dissoci- 

•  * 

ated  gases  and  related  phenomena’®’ with  application  to  rocket  motor  heat  trans- 

»  • 

*  fer.  The  results  of  these  studies  are  reported  elsewhere’^.  The  experimental  work 
carri^  out  thus  far  is  described  below. 

Because  of  the  dominant  influence  of  density  level  on  the  rate  of  chemical 
reactions  in  the  gas  phase,  it  was  decided  at  the  outset  to  carry  out  experfcients  over 
a  wide  range  of  pressures.  The  difficulties  of  operating  steady  flow  combustion  de¬ 
vices  of  reasonable  size  at  pressures  in  the  range  of  1000  psia  suggested  the  use  of 
a  constant  volume  transient  explosion  method  for  the  high  prelsure  studies.  However, 
at  the  lower  end  of  the  pressure  scale,  the  steady  flow,  flat  flame  subatmospheric 
burner  showed  great  promise.  Thus,  the  experimental  program  wa^ divided  into  two  dis¬ 
tinct  phases  which  we  will  arbitrarily  designate  hereafter  as  "superatmospheric"  and 
"subatmospheric".  In  the  superatmospheric  program  we  set  out  to  deve*lop  an  explosion 
technique  for  studying  transient  energy  transport  rates  in  reacting  gases  up  to  pres¬ 
sures  of  about  1000  psia.  In  the  subatmospheric  program  we  set  out  to  develop  tech¬ 
niques  for  studying  the  steady  state  heat  transfer  rate  to  <yoled  probes  immersed  in 

the  flame  gases  at  pressures  as  low  as  the  mm  Hg  range  (between  10"^  and  10"^  atm). 

0 

Both  phases  are  described  below,  and  some  of  the  experimental  results  obtained  to  date 
are  given.  It  is  felt  that  the  techniques  described  herein  have  yet  to  bear  the  larger 
portion  of  the  fruits  of  which^they  are  capable.  This  research  is  a  portion  of  a  con¬ 
tinuing  program  at  AeroChem  devoted  to  the  subject  of  energy  transfer  in  the  presence 
of  thermochemical  chang’e. 


II.  MEASUREMENT  OF  HEAT  TRANSFER  TO  SOLIDS  IiMMERSED  IN  FLAME  GASES 


AT  SUBATMOSPHERIC  PRESSURES 


1.  Nature  of  Problem 

The  work  of  Kilham’^’'''^’’^  on  radiation  cooled  circular  cylinders  in  cross 
flow’ indicated  that  anomalous  heat  transfer  effects  are  encountered  near  the  re- 
action  zone  of  flames  burning  at  atmospheric  pressure.  In  hydrogen-oxygen  flames  sta¬ 
bilized  on  a  Bunsen  type  burner,  for  example,  he  found  that  heat  transfer  rates  in  or 
just  above  the  inner  cone  were  higher  than  computed  values  based  on  "conventional" 
convection  correlations.  More  recently,  Cookson’*  and  Kilham  have  corroborated  these 
preliminary  findings  in  a  more  detailed  study  of  heat  transfer  from  flat  lean  hydrogen- 
air  flames  to  a  stagnation-point  heat-flow  meter;  also  at  atmospheric  pressure.  From 
this  work,  as  well  as  past  experience  with  the  immersion  pyrometry  of  flame  gases  in 
regions  of  active  chemical  reaction’’’’,  it  appears  that  this  enhancement  in  heat  trans¬ 
fer  rate  is  associated  in  some  way  with  the  presence  of  energetic  molecular  fragments 
(free  radicals)  in  the  gases  wetting  the  heat  transfer  surfaces.  Since  temperature 
and  heat  flux  measurements  are  often  found  to  be  sensitive  to  the  nature  of  the  sur¬ 
face  used  (e.g. ,  whether  it  is  a  bare  metal  or  coated  with  quartz),  the  postulated 
mechanism  of  diffusion  and  surface  catalyzed  chemical  reassociation  is  quite  plausible; 
although  the  effects  of  gas  phase  chemical  reaction  within  the  diffusion  boundary 
layer  cannot  be  ruled  out,  particularly  at  atmospheric  pressure. 

If  one  now  considers  low  pressure  flames,  and  in  particular  low  pressure 
flat  flames,  then  the  space  coordinates  of  the  combustion  process  are  considerably 
expanded  --  reaction  zone  thicknesses  are  measured  in  centimeters  rather  than  mic¬ 
rons’®’’^.  Thus,  while  the  dimensions  of  any  heat  transfer  probe  are  apt  to  be  large 
compared  to  the  reaction  zone  thickness  at  one  atmosphere,  in  a  subatmospheric  flame 
one  has  the  possibility  of  truly  probing  the  reaction  zone.  This  is  advantageous  in 
the  present  case  since  there  is  now  ample  reason  to  believe  that  free  radical  concen¬ 
trations  show  maxima  wi thin  the  reaction  zone^°.  Moreover,  these  maxima  are  far  in 
excess  of  the  equilibrium  free  radical  concentration  corresponding  to  complete,  adi- 
abatic  combustion^’.  For  simple  free  radicals,  the  dominant  gas  phase  removal  proc¬ 
esses  are  termolecular.  Thus,  one  has  the  further  advantage  that  the  nonequilibrium 
free  radical  concentrations  should  persist  further  downstream  of  the  primary  reaction 
zone  at  the  lower  pressures.  A  closely  related  phenomenon  is  the  reassociation  of 
free  radicals  in  the  vicinity  of  the  probe  interface.  Since  surface  catalysis  effects 
will  only  be  observed  under  conditions  of  long  radical  lifetime  compared  to  the 


residence  time  within  the  diffusion  boundary  layer,  there  are  a  priori  reasons  to  be¬ 
lieve  that  heat  transfer  specificity  effects  with  various  heat  transfer  surfaces  will 
become  more  and  more  marked  as  the  pressure  level  is  reduced.  This  trend  has  already 
been  reported  by  Cutting,  Fay,  Hogan  and  Moffatt^^  in  a  series  of  shock-tube  experi¬ 
ments  with  stoichiometric  oxy-hydrogen  mixtures.  • 

The  previous  considerations  suggested  the  construction  of  a  low  pressure 
flat  flame  burner  and  a  steady  state  heat  transfer  probe  which  would  enable  the  mech¬ 
anism  of  energy  transport  to  be  studied  at  various  points  within  the  flame.  Except 
for  minor  departures,  the  burner  design  was  patterned  after  similar  burners  which  have 
been  previously  operated  at  this  laboratory  in  connection  with  the  study  of  ion  pro¬ 
duction  in  flames^^.  However,  the  heat  transfer  probe  had  several  novel  aspects,  most 
of  which  are  described  in  the  following  section.  In  addition  to  heat  flux  measurements, 
some  work  has  been  done  on  gas  temperature  measurements  using  both  a  sodium  line  rever¬ 
sal  technigue,  and  a  compensated  hot  wire  method^'^’^^.  As  yet,  neither  of  these  meas¬ 
urements  has  been  made  on  a  routine  basis,  and  more  complete  results  of  this  work  must 
of  necessity  await  a  future  communication.  Clearly,  if  heat  flux  measurements  are 
ultimately  to  be  interpreted  in  terms  of  energy  transport  mechanisms,  determination 
of  the  heat  transport  coefficient  will  be  necessary.  For  this  purpose  measurements 
of  flame  gas  temperature  and  mass  velocity  will  be  a  necessary  (but  not  sufficient) 
adjunct  to  the  heat  flux  data.  Also  necessary,  but  more  distant,  would  be  an  indepen¬ 
dent  determination  of  the  free  radical  concentration  distribution  in  the  particular 
flames  under  study. 

2.  Heat  Transfer  Probe 


^A  right  circular  cylinder  oriented  transverse  to  the  flow  of  flame  gases 
can  form  the  basis  of  a  simple  calorimetric-type^^  heatgflux  probe,  as  shown  in  Fig.  1. 
In  this  design,  the  temperature  rise  Tj  -  T^  experienced  by  a  known  flow  rate  of  cool¬ 
ant  passed  through  the  tube  is  sufficient  to  determine  thi  average  heat  flux,  provided 
^that  edge  effects  in  the  flame  are  avoided  and  axial  heat  conduction  within  the  probe 
0is  rendered  negligible.  The  edge  effects  can  be  minimized  by  sensing  only  the  temper¬ 
ature  rise  over  a  length  L  of  tube  which  is  well  imbedded  within  the  flame  (i.e. ,  L  is 
considerably  smaller  than  the  burner  diameter).  If  axial  conduction  losses  wei^e  sig¬ 
nificant,  the  effects  could  always  be  experimentally  eliminated  from  the  data  by  ob¬ 
serving  the  temperature  rise  data  over  a  range  of  coolant,mass  flows,  thereby  enabling 

extrapolation  to  vanishing  temperature  ris^.  In  practice,  if  the  Peclet  number  Ul/a 

« 

based  ontaverage  coolant  velocity  U,  length  L,  and  coolant  thermal  diffusivity  a  is 
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sufficiently  large,  convective  energy  transport  will  dominate,  and  axial  conduction 

0 

at  the  thermocouple  locations  can  be  neglected  without  serious  error. 

^  A  further  design  considefation  has  to  do  with  how  representative  of  the  cool¬ 

ant  bulk  temperature  are  the  measurements  and  T2.  Under  conditions  of  high  heat 
flux  from  the  flame  gases  to  the  coolant,  the  radial  temperature  differences  in  the 
coolant  can  become  comparable  to  or  exceed  the  desired  axial  temperature  rise,  depend¬ 
ing  on  the  choice  of  coolant.  Calculations  showed  that  this  was  certainly  true  for 
gaseous  coolants  and  was  likely  to  be  true  for  ordinary  liquids,  such  as  water.  This 
means  that  a  small  error*in  the  in  the  axial  alignment  of  the  thermocouples  can  cause 
a  large  systematic  erjor*in  the  calculated  heat  flux.  Furthermore,  the  more  non-uniform 
are  the  radial  temperature  profiled  at  stations  1  and  2,  the  less  likely  are  the  center- 
line  temperature  difference  measurements  to  be  indicative  of  the  axial  rise  in  mean 
"mixed"  temperature.  For  these  reasons,  it  was  decided  to  operate  the  probe  with  a 
liquid  metal  coolan'f,  and  mercury  proved  to  be  a  convenient  choice.  With  this  choice, 
under  operating  conditions,  the  radial  temperature  differences  are  small  compared  to  the 
axial  temperature  rise  and,  furthermore,  a  reasonably  accuftite  estimate  of  the  outer 
surface  temperature  <ff  the  heat  transfer  probe  becomes  possible  from  a  knowledge  of  the 
observed  heat  flux  and  an  estimate  of  the  internal  (liquid  metal  side)  heat  transfer  co¬ 
efficient. 

The  overall  configuration  of  the  heat  transfer  probe  is  shown  schematically 
(plan  view)  in  Fig.  2.  Liquid  mercury  enters  the  probe  from  a  constant  head  flow  sys¬ 
tem  (to  be  described)  and  is  directed  through  tubing  into  a  crescent  'shaped  header 
straddling  the  burner.  The  mercury  then  flows  through  the  thin  walled  tubing  compris¬ 
ing  the  heat  transfer  surface  and  leaves  through  the  elbow  at  the  other  end  of  the  cres¬ 
cent  shaped  header.  Type  304  stainless  steel  was  used  for  all^parts  wetted  by  the 
liquid  mercury.  The  heat  transfer  tube  itself  was  fabricated  from  .125"  O.D. ,  .003" 
wall  thickness  seamless  tubing  supplied  by  J.  Bishop  and  Company  Platinum  1/forks,  Malvern 
Pennsylvania.  Fig.  3  shows  a  detailed  cross-section  of  the  probe  thermocouple  and  cool- 
ant  flow  configuration.  The  thermocouples  used  were  of  .040"  O.D.  stainless  steel 
sheathed  iron-constantan  and  were  supplied  by  Aero  Research  Company,  Chicago  7,  Illinois 
A  helix  of  stainless  steel  wire  was  used  to  maintain  co-axial  alignment  between  the 
thermocouples  and  the  heat  transfer  tube  itself.  The  packings  shown  in  black  are  of 
neoprene  and  are  either  standard  0-rings  or  standard  Conax  Company  (Buffalo,  New  York) 
packings.  In  earlier  versions  of  the  probe,  0.125"  O.D.  sheathed  thermocouples  had 
been  used  with  bare  junctions.  In  order  to  record  differential  temperature  rise,  how¬ 


ever,  these  jfinctions  had  to  be  coated  (electrically  insulated)  with  an  epoxy  type 


resin;  however,  this  proved  unsatisfactory  in  that  the  coating  life  in  the  heated  mer- 
cury  was  erratic.  Replacement  with  ungrounded  .040"  sheath  enclosed  thermocouples  eli¬ 
minated  the  necessity  of  this  coating  process  and  reduced  the  "out  of  service"  time  of 
the  probe  considerably.  For  most  of  the  work  reported  here,  the  burner  diameter  was 
10.1  cm  and  the  length  L  between  thermocouple  junctions  was  5. OS  cm. 

As  shown  in  Fig.  3,  heat  transfer  surfaces  could  be  interchanged  by  first 
loosening  the  appropriate  neoprene  and  packings.  This  was  considered  desirable  for 
studying  the  effects  of  surface  catalysis,  since  the  same  basic  probe  could  be  used 
for  heat  transfer  experiments  with  surfaces  of  widely  different  catalytic  activity. 

For  this  purpose  several  sections  of  the  thin  walled  stainless  tubing  were  electro- 
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plated  with  platinum  ,  and  several  sections  were  given  a  coat  of  silicone  (Dow  Corning 
No.  803).  Comparative  data  using  both  types  of  surface  were  then  taken  in  an  ethylene- 
air  flame  and  a  hydrogen-oxygen  flame  in  an  effort  to  observe  any  systematic  effects 
due  to  preferential  free  radical  recombination  on  the  platinum.  With  further  improve¬ 
ments  in  the  probe  technique,  we  hope  to  extend  the  limited  amount  of  data  presented 
herein  to  other  mixture  ratios  and  other  surfaces. 

The  arrangement  of  the  constant  head  mercury  flow  system  is  shown  in  Fig.  4. 
Basically,  mercury  flows  from  an  elevated  reservoir  (item  2)  through  the  probe  to  a 
Sigma-motor  "finger  pump"  (item  4)  which  continuously  replenishes  the  reservoir.  Early 
flow  rate  determinations  were  made  using  the  steady  stat^pressure  drop  across  a  re¬ 
striction  (item  5)  in  the  flow;  however,  this  method  was  abandoned  in  favor  of  a  simple 
volumetric  "batch"  method  (item  7)  owing  to  the  need  for  frequent  recalibration.  With 
this  coolarft  flow  system,  extreme  care  had  to  be  exercised  to  make  sure  that  small  air 
pockets  were  not  trapped  within  the  probe  itself  during  a  run. 

To  examine  experimentally  the  severity  of  axial  conduction  errors,  a  series 
of  bench  tests  was  run  at  atmospheric  pressure  using  a  shielded  Meker  burner.  Appar¬ 
ent  heat  flux  was  plotted  against  the  differential  emf  E2  -  E-j  over  a  range  of  coolant 
flow  rates.  In  a  typical  case,  more  than  a  two-fold  change  in  coolant  temperature  rise 
(56°C  to  about  122°C)  produced  no  systematic  change  in  apparent  heat  flux,  the  latter 
remaining  constant  to  within  +l^%.  The  Peclet  numbers  for  these  extremes  (based  on  the 
length  L)  were  550  and  726  respectively.  In  practice,  the  subatmospheric  flame  tests 
were  run  at  mercury  flow  rates  corresponding  to  E2  -  E,  ~  1  mv  or  about  20°C  temper¬ 
ature  difference^.  The  probe  data  recorded  included: 

^We  are  indebted  to  the  Chemical  Division  of  Engelhard  Industries,  Newark,  N.J.,  for 
^this  service.  •  • 

For  iron-constantan  in  the  temperature  range  investigated,  the  Leeds  and  Northrop 
Company  thermocouple  tables  give  1  mv  =  IS. 8  C,  to  within  +  1.51^. 
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1.  the  emf  Ei , 

2.  the  differential  emf  Ej  -  Ei ,  and 

3.  the  time  to  fill  a  burette  of  known  volume  (usually 
50  cc)  with  mercury. 

From  these  data,  the  average  heat  flux  d|"  was  computed  from  the  steady  state  energy 
balance  relation  • 


m,,  c  ,,  (at),, 
Hg  p,Hg  ^Hg 


7r  d  L 


where 


m  =  mass  flow  rate  (gm/sec)  of  mercury  through  probe 
^  (product  of  density  and  volume  flow  rate) 

c  =  average  specific  heat  of  mercury  [  (cal)  (gm)"U°C)"’] 

"  over  the  temperature  interval  AT 

(AT)„  =  axial  temperature  rise  of  mercury  (^C) 

Hg 

d  =  outer  diameter  of  heat  transfer  surface  (cm) 

L  =  axial  distance  between  thermocouple  junctions  (cm) 

For  this  purpose  the  thermal  and  physical  data  assembled  in  Appendix  1  were  used. 


3.  Subatmospheric  Burner 

The  burner  assembly  used  in  the  present  work,  shown  in  Fig.  5>  is  basicafty 
similar  to  those  used  for  earlier  flame  studies  at  this  laboratory^^.  The  burner  it¬ 
self  is  a  thin  walled, 10  cm  diameter  brass  cylinder  mounted  on  a  2.54  cm  steel  support 
and  supply  tube.  Within  the  brass  cylinder  are  a  series  of  five  stainless  steel  screens 
separated  by  spacers  for  uniformization  of  the  flow  of  pre-mixed  gases.  The  burner 
canTe  moved  vertically  by  means  of  a  threaded  section  of  the  support  tube  and  suit¬ 
able  gearing.  The  burner  support  conduit  comes  through  a  Teflon  packing  gland  in  the 
vacuum  jacket.  Since  the  burner  can  be  raised  and  lowered  while  in  operation,  one 
thereby  obtains  a  traverse  using  probes  which  are  fixed  v/ith  respect  to  the  burner 
housing.  The  subatmospheric  pressure  is  maintained  by  the  action  of  a  Stokes  mech¬ 
anical  vacuum  pump  of  nominal  pumping  capacity  360  cfm.  Four  large  openings,  at  right 
angles  to  one  another,  are  provided  in  the  housing  for  visibility  and  instrumentation 

9. 

lead-throughs.  The  housing  itself  was  fabricated  from  stainless  sheet  steel  and  pro- 
•  vided  with  four  slip  flanges  and  two  weld-on  flanges.  The  position  of  the  probe 
(mounted  by  means  one  of  the  slip  flanges)  with  respect  to  the  burner  screen  was 
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determined  from  the  reading  of  a  digital  counter  linked  to  the  gear  system  (l  count  = 
4.55  X  10"^  cm). 

Gas  flows  are  metered  by  individual  rotameters  (variable  area  flow  meters) 
together  with  a  total  gas  flow  rotameter.  The  mixture  ratio  is  computed  from  the  rea¬ 
dings  of  the  individual  rotameters;  however,  the  total  flow  is  not  the  sum  of  the  in¬ 
dividual  flows  because  of  the  presence  of  a  bleed  check  valve  introduced  between  the 
total  and  individual  flow  meters.  By  exhausting  varying  quantities  of  the  mixture  to 
the  atmosphere,  the  check  valve  enables  the  combustible  gas  composition  to  be  main¬ 
tained  nearly  constant  while  varying  either  the  total  flow  rate  or  pressure.  Thus, 
any  particular  flame  can  be  ignited  at  a  convenient  pressure,  and  the  pressure  level 
can  then  be  reduced  vi/hile  maintaining  the  mixture  composition  substantially  constant. 

Ignition  is  accomplished  by  passing  a  spark  between  an  electrode  and  the 
pilot  rim  in  the  presence  of  a  small  amount  of  pilot  gas,  the  latter  being  introduced 
through  a  separate  tube  of  small  bore.  The  resulting  diffusion  flame  at  the  burner 
rim  then  ignites  the  main  pre-mixed  gas  flow. 

4.  Experimental  Results 

To  check  out  the  burner-heat  flux  probe  system  as  a  whole,  the  first  series 
of  low  pressure  experiments  was  done  with  ethylene-air  flames  over  a  range  of  pressure 
and  mixture  ratio.  A  particular  flame  was  then  singled  out  in  order  to  study  the  re¬ 
peatability  of  the  measurements;  this  was  a  lean  ethylene-air  flame  (*  =  0.838)  at  a 
nominal  pressure  level  of  14  mm  Hg.  Fig.  6  shows  the  heat  flux  distributions  measured 
with  both  silicone  coated  and  platinized  probes  on  four  differenl;  days.  It  is  obser¬ 
ved  that; 

a.  peak  heat  fluxes  for  both  the  platinized  and  silicone  coated 
probes  are  nearly  the  same.  The  platinized  probe  peak  heat 
flux  is  slightly  larger  (about  %)  and  tends  to  occur  somev;hat 
upstream  of  that  for  the  silicone  coated  probe,  but  neither  of 
these  effects  is  particularly  definitive; 

b.  peak  heat  fluxes  for  both  probes  occurred  noticeably  down¬ 
stream  of  the  region  of  visible  luminousity.  In  the  case  of 
the  platinized  tube,  the  distance  downstream  amounted  to  about 
0.7  visible  flame  thicknesses  (measured  from  rear  of  zone). 

These  visible  flame  zone  thicknesses  were  about  2.5  times  the 
probe  diameter  (d  =0.318  cm); 
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c.  peak  i^eat  fluxes  were  about  0.75  (cal) (cm)"^ (sec)"^  at  this 
pressure  level.  The  heat  fluxes  decreased  with  distance  in 

the  burned  gases  at  about  the  same  rate  for  both  probe  sur-  * 
faces.  The  major  portion  of  this  decrease  is  attributed  to 
radial  heat  loss,  perhaps  augmented  by  recirculation  flows 
within  the  burner  housing  downstream  of  the  luminous  zone; 

d.  for  any  given  run,  experimental  heat  fluxes  at  a  prescribed 
distance  from  the  burner  could  be  reproduced  to  within  about 
+3^.  On  different  days,  however,  the  reproducibility  is  of 
the  order  of  only  +15JS  in  the  regions  of  small  (rapidly  ris¬ 
ing)  heat  flux,  upstream  of  and  within  the  luminous  zone. 

The  major  cause  of  this  loss  in  precision  is  felt  to  be  var¬ 
iation  in  the  position  of  the  flame  with  reject  to  the  burner 
rim.  This  will  be  discussed  further  in  connection  with  re¬ 
lated  experiments  in  the  hydrogen/oxygen  system. 

Since  free  radical  concentrations  in  hydrocarbon-air  flames  are  believed  to  be  quite 
small^"^,  the  ethylene-air  results  summarized  above  are  not  unexpected.  Because  of 
the  complexity ^f  the  gas  mixture  enveloping  the  probe  surfaces  (composed  of  CO,  COjj 
O2,  Mjj  HjO  and  smaller  amounts  of  OH  and  O),  theoretical  interpretation  of  the 
ethylene-air  flame  results  is  difficult.  For  these  and  other  reasons,  the  experimen¬ 
tal  effort  was  redirected  to  the  hydrogen-oxygen  flames,  in  which  there  is,  at  present, 
considerable  theoretical,  as  well  as*practical,  interest. 

•  ♦  In  the  case  of  the  hydrogen-oxygen  flame,  the  primary  reaction  zone  is  not 

visible.  Furthermore,  as  suggested  by  the  ethylene-air  flame  data  discussed  above, 
the  distance  from  the  probe  to  the  burner  screen  is  not  a  satisfactory  measure  of  the 
probe  location  within  the  flame  zone  since  the  flame  "stand  off  distance"  may  vary 
from  day  to  day.  One  thus  requires  a  position  measurement  relative  to  the  reaction 
zone  itself'.  Two  general  methods  come  to  mind,  viz,  an  immersion  method,  such  as 
the  use  of  an  auxiliary  thermocouple  fixed  with  respect  to  the  probe,  or  an  external 
radiation  detector  fixed  with  respect  to  the  probe.  Even  if  the  reaction  zone  were 
visible,  the  human  eye  is  inadequate  for  accurately  determining  the  flame  position. 

For  the  hydrogen-oxygen  flame,  the  human  eye  is  out  of  the  question.  It  was  therefore 

'it  will  probably  be  desirable  to  control  the  burner  rim  temperature,  as  v/ell,  in 
future  experiments. 

•  » 

M, 
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decided  to  employ  a  large  diameter  radiation  cooled  Pt/Pt-13^  Rh  thermocouple  wire 
(ceramic^  coated  to  a  nominal  diameter  of  2.5  mm)  stretched  parallel  to  th<3  probe 
surface  with  the  junction  located  at  the  mid  point.  By  shifting  the  observed  heat 
flux  data  in  the  streamwise  direction  until  the  "temperature"  profiles  obtained  with 
the  thermocouple  were  aligned,  the  relative  output  of  the  thermocouple  could  then  be 
used  to  ascertain  the  probe  position  within  the  flame  zone.  This  procedure  has  the 
advantage  that  both  the  heat  flux  and  "temperature"  can  be  recorded  simultaneously 
with  no  particular  difficulty.  * 

Fig.  7  shows  a  typical  pair  of  such  profiles  obtained  with  a  platinized 
probe  in  a  rich  ('^  =  1.43)  oxy-hydrogen  flame  at  5  mm  Hg  pressure.  A  striking  feature 
evident  here,  as  in  all  of  the  5  mm  Hg  oxy-hydrogen  data  thus  far  exaiftned,  is  the 
rapid  spatial  decay  of  heat  flux  in  the  downstream  portions  of  the  flame;  i . e. ,  the 
heat  flux  decays  at  a  rate  comparable  to  the  rate  of  increase  in  the  upstream  regions 
of  the  flame.  Also  of  interest  is  the  fact  that  the  point  of  maximum  apparent  tem¬ 
perature  occurs  noticeably  downstream  of  the  point  of  peak  heat  flux.  This  was  also 
observed  for  silicone  coated  surfaces  in  the  same  flame,  and  is  suggestive  of  energy 
transfer  mechanisms  other  than  ordinary  heat  convection  within  the  reaction  zone.  A 
total  of  6  traverses  were  made  in  this  flame  (3  using  platinum  probe  surfaces,  3 
using  silicone  coated  surfaces)  in  the  hope  of  being  able  to  study  the  differences 
between  the  heat  flux  to  each  of  these  surfaces  at  a  function  of  position.  Surpris¬ 
ingly  enough  (in  the  light  of  the  ethylene-air  data)  the  repeatability  proved  so  poor 
that  definitive  conclusions  could  not  be  drawn.  Thus,  while  the  maximum  observed 
value  of  the  peak  heat  flux  with  a  platinum  coated  probe  exceeded  the  maximum  value 
of|^:he  peak  heat  flux  with  a  silicone  coated  surface  by  as  much  as  17  percent,  the 
minimum  recorded  peak  heat  flux  with  a  platinum  probe  was  less  than  the  minimum  recor 
ded  silicone  coated  probe  peak  heat  flux  by  about  this  same  amount.  As  of  the  moment 
the  cause  of  this  lack  of  repeatability  in  the  oxy-hydrogen  system  can  only  be  spec¬ 
ulated  upon  (see  part  II-5).  Future  research  will  therefore  be  directed  at  this  most 
immediate  obstacle  to  the  elucidation  of  gas/solid  ener^  transfer  mechanisms  in  chem 
ically  reacting  mixtures  at  subatmo spheric  pressures. 

*^using  Astroceram  "A"  high-temperature  cement  (American  Thermo  catalytic  Corp.,  New 
York,  manufacturers) 

^particularly  in  the  case  of  the  platinum  probe  data,  v^ere  the  peak  heat  fluxes  on 
3  different  day^varied  by  as  much  as  about  +15^.  The  value  shown  in  Fig.  7  was 
the  lowest  in  absolute  value. 
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5.  Summary  of  Subatmospheric  Phase 

A  novel  method  of  measuring  steady-state  heat  fluxes  to  strongly  cooled 
catalytic  and  non-cata lytic  solids  within  subatmospheric  flames  (in  the  mm  Hg  range) 
has  been  described  and  used  to  obtain  preliminary  heat  flux  data  as  a  function  of 
position  in  flat  ethylene/air  and  hydrogen/oxygen  flames  stabilized  on  a  4  inch  dia¬ 
meter  burner.  The  present  probe  lends  itself  to  further  miniaturization,  (The  Bishop 
Company  supplies  thin  walled  stainless  tubing  smaller  by  a  factor  of  two  than  that 
used  here)  and  can  be  operated  for  long  periods  of  time  over  a  wide  range  of  total 
pressures. 

In  the  case  of  a  lean  -  .84)  ethylene/air  flame  at  14  mm  Hg  total  pres¬ 
sure  peak  heat  fluxes  to  both  platinized  and  silicone  coated  probe  surfaces  occurred 
noticeably  downstream  of  the  region  of  visible  luminousity.  However,  no  appreciable 
difference  in  peak  heat  flux  to  the  different  surfaces  was  observed  for  this  hydro¬ 
carbon/air  flame.  Similar  measurements  in  a  rich  oxy-hydrogen  flame  at  5  mm  total 
pressure  suggest  that  the  region  of  peak  heat  flux  to  both  platinized  silicone  coated 
probe  surfaces  occurs  upstream  of  the  region  of  maximum  gas  temperature.  Here,  how¬ 
ever,  the  precision  of  the  data  was  not  adequate  to  quantitatively  study  the  depen¬ 
dence  of  the  difference  in  heat  fluxes  as  a  function  of  position.  The  wide  scatter 
of  the  data  for  the  platinum  coated  probe  suggests  that  the  reproducibility  of  the 

t 

platinum  surface  itself  may  have  to  be  controlled  more  carefully  than  was  done  here  . 
Future  research  would  be  directed  at  this  obstacle  and  have  among  its  foremost  goals 
the  acquisition  of  data  on  the  dependence  of  the  difference  in  heat  flux  on  both 
position  and  total  pressure  within  chemically  reacting  hydrogen-oxygen  mixtures. 
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Platinized  probe  surfaces  were  degreased  with  trichloroethylene  prior  to  each  run. 
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III.  TRANSIEOT  HEAT  TRANSFER  RATES  FROM  THE  HIGH  PRESSURE  AND  TEMPERATURE 


PRODUCTS  OF  hydrogen-oxygen  EXPLOSIONS  IN  A  VESSEL  OF  CONSTANT  VOLUME 

1.  Introduction 

Since  the  experimental  complexity,  and  therefore  the  cost  of  actual  liquid 
propellant  rocket  motor  heat  transfer  tests,  is  so  large  one  is  led  to  inquire  whether 
something  can  be  learned  about  the  energy  transfer  mechanisms  involved  without  having 
to  operate  actual  rocket  motors  at  pressures  up  to  1000  psia.  Preferably,  the  experi¬ 
mental  method  sought  should  have  the  attributes  of  low  cost  and  simplicity;  i.e. ,  it 
should  enable  rapid  testing  in  the  laboratory  and  utilize  small  amounts  of  propellant 
without  the  need  for  extensive  facilities  or  equipment.  A  similar  question  has  re¬ 
cently  been  answered  in  the  affirmative  by  Griffin,  Turner  and  Angeloff^®  with  regard 
to  the  determination  of  the  specific  impulse  of  new  propellants.  These  investigators 
have  shown  that  constant  volume  combustion  in  a  ballistic  bomb  can  provide  accurate 
data  from  which  the  specific  impulse  can  be  obtained.  We  propose  here  to  show  that 
certain  important  molecular  transport  coefficients  (viz,  the  thermal  diffusivity  and 
radiative  emissivity)  can  be  profitably  investigated  in  a  constant  volume  bomb.  This 
type  of  information  would  not  only  comprise  a  useful  supplement  to  experimental  data 
obtained  in  flow  systems  using  the  same  propellants  but  also  could  be  incorporated 
into  a  more  comprehensive  theoretical  analysis  of  energy  transfer  processes  in  chem¬ 
ical  reactors  in  general  (and  rocket  thrust  chambers  in  particular). 

In  what  follows,  we  will  discuss  certain  theoretical  and  experimental  as¬ 
pects  of  the  use  of  constant  volume  explosions  to  study  energy  transport  in  chemi¬ 
cally  reacting  systems. 

2.  Theory 

Explosions  of  premixed  gases  in  closed  vessels  have  been  studied  for  many  ^ 
years^^’^°;  however,  the  emphasis  has  been  upon  the  reaction  time  and  relation  be¬ 
tween  the  peak  pressure  and  the  thermodynamic  properties  of  the  reactants  and  pro¬ 
ducts.  The  recent  work  of  Medvedev^’  (1958)  has  once  again  made  it  clear  that  spher¬ 
ical  bomb  explosion  studies  can  yield  very  useful  results  pertaining  to  the  equilibrium 
thermodynamic  properties  of  gases  at  high  temperatures  and  pressures.  It  appears  to 
the  writer  that  the  potentialities  of  this  technique  for  determining  transport  pro¬ 
perties  of  high  temperature  and  pressure  gases  remain  largely  unexploited. 
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The  essentials  of  the  explosion  method  can  be  described  as  follows.  A  pre¬ 
mixed  combustible  gas  mixture  is  admitted  into  a  pressure  vessel  cavity  and  ignited 
at  the  center  of  this  "bomb"  by  either  a  spark  or  an  electrically  heated  wire.  The 
chemical  energy  release  following  the  ignition  process  causes  the  pressure  in  the 
vessel  to  rise  rapidly  to  some  peak  value  p^  in  a  time  t^  at  which  the  reactants  may 
be  said  to  be  consumed.  At  this  point  the  hot  explosion  products,  which  are  now  in 
contact  with  the  cooler  walls  of  the  bomb,  begin  to  lose  appreciable  quantities  of 
heat,  principally  by  conduction  and  radiation.  Since  the  internal  energy  of  the  en¬ 
closed  gases  is  related  to  the  pressure,  the  cooling  period  is  accompanied  by  a  gradual 
pressure  decay  (see  Fig.  8)  which  is  easily  measured.  Many  good  records  of  this  type 
may  be  found  in  the  literature  and,  as  will  be  seen,  can  be  interpreted  from  the  point 
of  view  of  energy  transfer  in  reacting  systems.  Unfortunately,  however,  some  of  the 
best  records  have  been  obtained  in  cylindrical  bombs  of  small  aspect  ratio  with  non¬ 
central  ignition.  In  the  past,  more  interest  has  been  displayed  in  the  pressure  rise 
portion  of  the  record  than  in  the  detailed  behavior  of  the  "cooling  curves".  For  the 
present  purpose,  this  emphasis  will  be  reversed. 

a.  Relation  Between  Average  Heat  Transfer  Rate  and  Rate  of  Pressure  Decay. 

To  appreciate  the  relation  between  the  pressure  history  and  the  amount  of  heat  loss, 
consider  a  mass  of  gas  of  volume  V  with  local  internal  energy  per  unit  mass  e  and 
local  density  ^.  Then  at  any  instant  t  of  time,  the  energy  content  of  the  gas  within 
the  vessel  will  be 

E(t)=JpedV  (2) 

If  the  specific  internal  energy  e  is  considered  to  be  a  function  of  p  and  T  alone, 
then  the  time  rate  of  change  of  E  can  be  written 


(3) 


Since  the  pressure  will  be  nearly  uniform  throughout  the  volume,  the  perfect  gas  law 
implies  that 


^  3 

dt 


(4) 
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Expanding  the  integrand,  it  is  observed  that,  under  th5  following  conditions^ 


'll  =  constant  =  'll 


(“|)  =  constant  =  (||) 

P  P 


~  _/8e\  „/9e\ 

e  -  T(g^)  -  P(9p) 

P  T 


(5a) 

(5b) 

(5c) 


dE/dt  and  dp/dt  will  be  linearly  related  by  the  simple  law 


dt  'r  ^  dt 

P  u 


(6) 


even  when  the  internal  energy  e  and  hence  the  temperature  T  vary  from  point  to  point 
within  the  vessel.  The  average  heat  transfer  rate  q"  per  unit  area  is  then  given  by 
(1/S) (dE/dt)  or 


q  gigj; 

P  u 


(7) 


where,  for  a  sphere, V/S  =  R/3.  The  total  heat  removed  up  to  time  t  (measured  from 
the  moment  of  peak  pressure)  is: 


Q(t)  =S 


r(t)dt=v(|f)  (f)[p  -p(t)] 

0  p  u 


(8) 


and,  in  particular,  Q(oo)  will  be  given  by  Eq. (8)  when  p(t)  is  evaluated  at  t  =  oo*. 


t"  9  G 

for  a  perfect  gas  e  =  e(T)  so  that  (^)  E  0.  Condition  (5c)  is  then  satisfied  for 

ap  j 

a  perfect  gas  with  constant  c  (specific  heat),  if  e(o)  =0. 

l^hen  water  vapor  condensation  occurs,  the  observed  pressure  p(oo)  must  be  corrected 
upwards  to  give  the  pressure  that  would  have  existed  had  no  condensation  taken  place. 
The  influence  of  condensation  is  discussed  further  in  Appendix  2. 
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b.  Graphical  Determination  of  Relative  Heat  Loss  Q  Using  Thermodynamic 
Charts  (Energy-Density  Diagram).  An  alternative  procedure  for  determining  Q  (t)  seems 
more  reasonable  when  applied  to  gas  mixtures  for  which  the  conditions  (5a) (5b) (5c) 
are  not  closely  satisfied.  Thus,  instead  of  taking  Q(t)  ~  ^p(t),  recourse  will  be 
made  to  a  graphical  representation  of  the  caloric  equation  of  state  e(p;p)  of  the 
products  of  reaction.  It  is  recalled  that  the  explosion  process  itself,  as  well  as 
the  subsequent  cooling  process, takejplace  in  a  constant  volume  with  no  loss  of  mass 
through  the  vessel  walls.  Therefore,  if  the  actual  non-uniform  gas  mixture  which  ex¬ 
ists  at  each  time  t  is  replaced  by  a  fictitious  uniform  gas  in  thermochemical  equili¬ 
brium  at  the  same  pressure  and  density,  the  cooling  process  for  this  fictitious  gas 
may  be  considered  to  be  represented  by  a  vertical  line  (line  of  constant  density)  on 
the  e-p  (see  Fig.  9)  plane.  Then  by  matching  the  observed  pressures  with  the  corres¬ 
ponding  isobars,  one  can  plot  the  normalized  heat  loss 

_Ae(t) 

^  “Ae(c») 

where  Ae(oo)  can  be  evaluated  from  the  experimental  pressures  p  and  p(oo)  or  obtain 

t  ™ 

ed  from  an  a  priori  knowledge  of  the  final  gas  temperature  . 

An  accurate  internal  energy-density  diagram  for  the  equilibrium  products 
of  stoichiometric  hydrogen-oxygen  explosions  is  shown  in  Fig.  10.  This  chart  was 
constructed  using  an  IBM  650  to  obtain  numerical  solutions  to  the  equations  of  mass 
action  and  overall  stoichiometry.  All  input  thermochemical  data  and  reference  ener¬ 
gies  are  those  of  Ref.  32. 

c.  Relation  Between  Relative  Heat  Loss  Q  and  Molecular  Transport  Proper- 
ties  of  the  Constituent  Gases.  Since  the  relative  heat  l(^s  Q  is  not  an  end  in  it- 
self,  we  must  inquire  further  as  to  the  relation  between  Q  (t)  and  the  energy  trans¬ 
port  properties  of  the  confined  high  temperature-high  pressure  gases. 


^If  the  heat  capacity  of  the  bomb  is  sufficiently  large,  the  final  gas  temperature 
(at  t  =  co)  will  be  very  nearly  T^,  the  initial  bomb  wall  temperature. 


To  illustrate  one  method  by  which  this  problem  might  be  approached,  consider 

f 

for  the  moment  that  all  of  the  heat  loss  occurs  by  conduction  to  the  cold  walls  and 
imagine  that  the  sphere  of  hot  gas  existing  at  the  moment  of  peak  pressure  behaves 
like  a  solid  sphere  of  constant  thermal  diffusivity  a,  "quenched"  by  sudden  contact 
with  cold  walls  at  radius  R.  Then  Q  would  increase  with  time  in  accord  with  the  well 
known  solution^^  to  the  Fourier  transient  heat  conduction  equation 

00 

(t:)  ~  ^  ~  ^  exp(-n2  %)  (ll) 

n=l 


where 

T  =  at/R^  =  nondimensional  time 
Q  (t)  =  normalized  integrated  heat  loss 
a  =  thermal  diffusivity 
R  =  radius  of  spherical  bomb 

Values  of  this  function  are  given  in  Table  I  and  plotted  in  Fig.  11.  Observe  that 
for  small  values  of  x  we  have  the  asymptotic  behavior 

Vt  (12) 

If  data  were  available  at  sufficiently  small  values  of  t  this  result  could  be  used 
to  obtain  experimental  values  of  the  effective  diffusivity  a  by  plotting  (7r/36)  (Q 
as  a  function  of  t/R^  and  measuring  the  slope  of  the  best  straight  line  through  the 
data. 

To  illustrate  this  procedure,  let  us  re-examine  the  conductive  heat  loss 
data  reported  by  David^'^  (1920)  for  a  \%  coal  gas-air  explosion  in  a  short  cylinder 
with  d  =  30  cm  and  l/d  =  1. 

In  constructing  Table  II,  we  have  corrected  David's  data  for  heat  losses 

before  the  time  of  peak  pressure  (i .e. ,  0.05  sec)  by^ subtracting»from  the  apparent 

value  of  Q  the  value  of  Q  at  t  =  t  .  Fig.  12  shows  a  plot  of  -^(Q  )^  against  t/R^. 

am  Jo 

The  inferred  thermal  diffusivity  in  this  case  is  about  4.8  cm^/sec. 


t 


Radiation  loss  will  be  discussed  in  a  subsequent  section. 


TABLE  I 


TRAMSIEOT  DIFFUSION  IN  A  SOLID  SPHE.RE^^ 


T 

1  -  Q 

* 

Q 

0.000 

1.000 

0.000 

0.00000 

0.001 

0.896 

0.104 

0.03162 

0.002 

0.355 

0.145 

0.04472 

0.003 

0.824 

0.176 

0.05477 

0.004 

0.798 

0.202 

0.06325 

0.005 

0.776 

0.224 

0.07071 

0.006 

0.756 

0.244 

0.07746 

0.007 

0.739 

0.261 

0.03367 

0.008 

0.721 

0.279 

0.08947^ 

0.009 

0.706 

0.294 

0.09487 

0.01 

0.691 

0.309 

0.1000 

0.02 

0.581 

0.419 

0.1414 

0.03 

0.504 

0.496 

0.1732 

0.04 

0.443 

0.557 

0.2000 

0.05 

0.393 

0.607 

0.2236 

0.06 

0.351 

0.649 

0.2449 

0.07 

0.314 

0.686 

0.2646 

0.08 

0.233 

0.717 

0.2828 

0.09 

0.254 

0.746 

0.3000 

0.10 

0.230 

0.770 

0.3162 

0.15 

0.139 

0.861 

0.3873 

0.20 

0.084 

0.916 

0.4472 

0.25 

0.052 

0.948 

0.5000 
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TABLE  11 


FRACTIONAL  HEAT  LOSS  BY  CONDUCTION;  [DAVID^^  (1920)' 
15^  Coal  Gas-Air  in  30  cm  L/d  =  1  Cylinder 


t  - 1 

0 

(sec) 

«a 

t  -  t 

m 

(sec) 

Q*  -  .051 

a 

0.05 

0.051 

0.00 

0.000 

0.10 

0. 144 

0.05 

0.093 

0.15 

0.203 

0.10 

0.152 

0.20 

0.245 

0.15 

0.194 

0.25 

0.277 

0.20 

0.226 

0.30 

0.300 

0.25 

0.249 

0.40 

0.344 

0.35 

0.293 

0.50 

0.376 

0.45 

0.325 

Now,  in  practice,  the  attainable  accuracy  of  a  method  of  this  type  appears 
to  be  limited  by  four  considerations: 

1.  the  pressure  and  temperature  dependence  of  the  thermal  diffusivity 
of  gases, 

2.  the  compressibility  of  the  gas  mixture  within  the  fixed  enclosure, 

3.  possible  enhancement  of  energy  transfer  due  to  gravity  induced 
natural  convection  currents  within  the  vessel,  and 

4.  energy  loss  by  other  mechanisms  (e.g. ,  radiation) 

The  first  two  of  these  factors  will  be  briefly  discussed  in  the  remainder  of  this 
section.  The  remaining  factors  will  be  dealt  with  in  separate  sections. 


Since  the  thermal  diffusivity  a  of  a  gas  depends  almost  linearly  upon  the 
inverse  of  the  pressure  level,  the  previous  procedure  can  only  be  useful  when  Ap/p^ 
is  small;  i .e. ,  in  effect  the  diffusivity  of  the  medium  increases  during  the  experi¬ 
ment.  Thus,  the  method  outlined  above  would  inherently  give  the  diffusivity  at  the 
pressure  p^  only  by  extrapolation  to  zero  time  (t  O).  This  is  not  a  serious  lim¬ 
itation  in  itself  since  p  can  be  varied  at  will  by  varying  the  initial  pressure 

m 

(before  explosion). 


Because  the  thermal  diffusivity  is  also  temperature  dependent,  a  method 
of  data  reduction  of  the  type  described  earlier  also  yields  an  "effective"  constant 
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property  diffusivity  which  after  extrapolation  to  zero  time  is  some  average  between 

srid  a(p  »T  ).  If  the  functional  form  of  the  temperature  dependence  of  a  were 
known,  then  solutions  to  the  transient  heat  conduction  equation  for  variable  diffu¬ 
sivity  could  be  coupled  with  the  experimental  data  to  yield  a(p  ,T  ).  For  instance, 

m  m 

Hansen's  one  dimensional  solutions^^*^®  to  the  heat  conduction  equation  could  be  car¬ 
ried  over  to  the  case  of  spherical  symmetry  with  no  difficulty.  An  interesting  ques¬ 
tion  which  arises  when  this  is  done,  however,  has  to  do  with  effects  of  thermally  in¬ 
duced  volume  changes  on  the  behavior  of  a  completely  enclosed  compressible  gas. 
Qualitatively  speaking,  when  the  exterior  regions  of  gas  are  suddenly  cooled,  a  con¬ 
traction  occurs  which  must  be  compensated  for  by  a  radial  expansion  of  the  central 
regions  of  hotter  gases.  The  result  of  this  expansion  process  is  to  increase  the 
temperature  gradient  at  the  wall,  and  hence  increase  the  rate  of  heat  loss  at  each 
time  t.  To  the  writer's  knowledge,  the  only  treatment  of  this  problem  is  the  linear¬ 
ized  treatment  of  Rayleigh^"^  (1899).  Rayleigh  clearly  stated  that,  even  in  the  limit 
of  infinitesimal  distrubances:  "llhen  one  part  of  the  gas  expands  in  consequence  of 
the  reception  of  heat  by  radiation  or  conduction,  it  compresses  the  remaining  parts 
and  these  in  their  turn  become  heated  in  accordance  with  the  law  of  gases.  To  take 
into  account  this  effect,  a  special  investigation  is  necessary".  His  solution  to  the 
linearized  problem  gives  Q  as  a  function  of  a  new  nondimensional  time  t'  related  to 
T  by  the  factor  Tr^/y,  where  y  =  c  /c  =  1.41  .  Wien  y  =  1,  the  solution  to  Rayleigh's 
problem  reduces  to  the  Fourier  solution  for  a  spherical  solid.  It  is  observed  that 
the  use  of  the  Fourier  solution  [Eq.(ll)]  when  in  fact  y  =  1.41  causes  the  inferred 
thermal  diffusivity  to  be  systematically  overestimated  by  about  40^.  This  points  up 
the  need  for  more  realistic  solutions  to  the  transient  (nonlinear)  heat  conduction 
problem  than  are  currently  available  in  the  literature.  As  is  apparent,  the  accuracy 
of  the  property  value  inference  depends  upon  the  accuracy  of  the  available  solutions. 
Both  the  Fourier  and  Rayleigh  solutions  are  strictly  applicable  to  the  case  of  con¬ 
stant  thermal  diffusivity.  It  can  be  verified  from  the  foregoing  analysis  of  David's 
data,  however,  that  most  of  the  available  points  fall  in  a  region  of  time  which  not 
only  exceeds  the  range  of  validity  of  the  asymptotic  relation  (l2)  but  also  exceeds 


'unfortunately,  Rayleigh's  solution  does  not  have  the  self-simi larity^property  su¬ 
ggested  by  the  inclusion  of  y  in  the  time  variable.  The  values  of  Q  versus  t  giv¬ 
en  in  Table  III  are  thus  valid  only  for  y  =  1.41  and  cannot  be  obtained  from 
Fig.  11  by  merely  re-interpreting  the  independent  variable. 


TABLE  in 


RAYLEIGH'S  SOLUTION^?  (l899);  ’TRANSIENT  HEAT  CONDUCTION  IN  A 
SPHERICAL  MASS  OF  AIR  (Y  =  1.4l)  CONFINED  BY  WALLS  AT  CONSTANT  TEMPERATURE 


X' 

1  -  Q 

* 

Q 

r 

yr 

0.000 

• 

1.000 

0.0000 

0.000000 

0.00000 

0.05 

0.7037 

0.2963 

0.007143 

0.08452 

0.10  , 

0.6037 

0.3963 

0.01429 

0.1195 

0.184 

0.5000 

0.5000 

0.0263 

0.162 

0.20 

0.4811 

0.5189 

0.02857 

0.1690 

0.30 

0.4002 

0.5998 

0.04286 

0.2070 

0.40 

0.3401 

0.6599 

0.05715 

0.2391 

0.50 

0.2926 

0.7074 

0.07143 

0.2673 

0.60 

0.2538 

0.7462 

0.08572 

0.2928 

0.70 

0.2215 

0.7785 

0.1000 

0.3162 

0.80 

0.1940 

0.8060 

0.1143 

0.3381 

0.90 

0.1705 

0.8295 

0.1286 

0.3586 

1.00 

0. 1502 

0.8498 

0.1429 

0.3780 

1.50 

0.0809 

0.9191 

0.2143 

0.4629 

2.00 

0.0441 

0.9559 

0.2857 

0.5345 

by  interpolation 


the  range  of  values  in  which  near  constancy  of  the  thermal  diffusivity  is  a  good 
approximation.  Since  the  thermal  diffusivity  a  appears  in  these  transient  conduction 
solutions  as  a  multiplier  of  the  time  variable  t,  the  effect  of  variations  in  a  is  to 
"stretch"  the  time  scale  of  the  process.  Using  this  notion,  the  full  Rayleigh  solution 
can  be  utilized  to  obtain  the  diffusivity  as  follows. 

Let  us  assume  that  the  effective  thermal  diffusivity  is  a  function  only  of 
the  pressure,  p,  (which  is  nearly  uniform  throughout  the  vessel  at  any  time)  and  of 
some  average  ("film")  temperature  T^.  Since  both  p  =  p(t)  and  T^  =  T^(t),  the  average 
diffusivity  can  be  regarded  as  some  function  of  time.  This  being  the  case,  inspection 
of  the  transient  heat  conduction  equation  reveals  that  the  constant  property  solution 
will  apply  to  this  situation  by  merely  re-interpreting  the  variable  t  as 


X 


a(t)  dt 


(13) 
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For  example,  as  a  reasonable  first  approximation,  consider  the  "hard-rsphere"  kinetic 
theory  dependence 

a  ~  p"’  (14) 

Let  us  further  assume  that  the  "film"  temperature  is  not  very  different 
from  the  arithmetic  mean  of  the  average  gas  temperature  T  in  the  vessel  (as  determined 
from  the  pressure)  and  the  wall  temperature  T  (nearly  constant),  i.e. 

W  ■ 

T^=i(T+Tj  (15) 

In  the  integrand  of  Eq.(l3)  we  then  have 

a(t)  =a^  .  (I,/I„)’/Mpyp)  *  (16) 

where  T^(t)  is  given  by  Eq.(l5).  If  we  define  a  new  "stretched"  time  tj^  by 

=  (*)>A  J  ‘  [(p/p„) .  (Tyi„)]’^'  {p/pj-  dt  (17) 

0 

then  T  =  (a  /R^)t,.  .  Since  the  integrand  of  Eq.(l7)  is  known  if  the  pressure  record 
fn  A 

p(t)  is  known,  the  calculation  of  'tj^(t)  presents  no  difficulty.  The  determination  of 
then  proceeds  as  follows.  From  a  knowledge  of  the  relative  conductive  heat  loss 
up  to  time  t,  the  corresponding  value  of  t  is  obtained  from  the  Rayleigh  solution. 
Values  of  t  obtained  in  this  manner  could  then  be  plotted  against  If  the  pos¬ 

tulated  heat  loss  mechanism  is  correct,  a  straight  line  of  slope  would  be  obtained. 

d.  Possible  Buoyancy  Effects;  Transient  Natural  Convection.  Another  impor¬ 
tant  question  in  the  use  of  the  constant  volume  explosion  technique  is:  How  much  time 
is  available  before  natural  convection  effects  would  significantly  disturb  heat  trans¬ 
fer  processes  occurring  within  the  vessel?  This  question  is  difficult  to  answer 
a  priori,  and  only  a  limited  amount  of  available  explosion  data  bear  upon  it.  No  use¬ 
ful  generalizations  have  appeared  in  which  the  role  of  the  important  parameters  have 
been  displayed. 

Lewis  and  Von  Elbe^®’^^  have  observed  that  under  certain  conditions  a  break 
appears  in  the  later  history  of  the  cooling  curve  at  which  point  the  heat  transfer  rate 
to  the  vessel  walls  is  suddenly  increased.  This  break  is  most  noticeable  if  the 
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explosion  products  have  a  large  thermal  conductivity  and  low  density;  substitution  of 
heavy  diluent  gases  (e.g. ,  argon)  for  light  diluents  (e.g. ,  helium)  makes  the  effect 
more  gradual.  Lewis  and  Von  Elbe  interpret  this  phenomenon  as  constituting  further 
evidence  that  a  radial  temperature  gradient  exists  in  the  gas  at  the  moment  of  peak 
pressure.  It  is  the  convective  rise  of  this  hot  core  of  gas  which  presumably  increase 
the  heat  transfer  rate.  Actually,  even  if  the  gas  temperature  at  the  moment  of  peak 
pressure  were  uniform,  it  is  possible  to  account  for  this  phenomenon  by  merely  invoking 
the  temperature  non-uniformities  brought  about  by  sudden  contact  with  the  cold  walls. 
Indeed,  it  will  be  seen  to  be  possible  to  make  an  estimate  of  when  this  convective 
augmentation  should  occur,  according  to  this  latter  model. 

In  his  study  of  steady  state  free  convection  in  horizontal  concentric  cyl¬ 
inders,  Kraussold'^°  found  that  there  was  a  threshhold  value  of  the  product 

Gr  •  Pr^  =  (b^gpAT/v^) (v/a)  (l8) 

beyond  which  the  effect  of  natural  convection  became  noticeable,  i .e. ,  no  appreciable 

bouyancy  effects  were  observed  until  Gr  •  Pn  >  10^.  Since  no  transient  data  applic- 

A 

able  to  the  sphere  is  available,  the  following  analogy  seems  reasonable  for  order  of 
magnitude  purposes. 

Imagine  the  transient  sphere  problem  to  be  equivalent  to  a  series  of  steady 
state  concentric  sphere  experiments  with  slowly  increasing  clearance.  Physically,  the 
clearance  over  which  the  large  temperature  change  occurs  becomes  the  thickness  of  the 
ther|?al  layer  which  propagates  into  the  hot  gas.  At  the  time  the  Grashof  number-Prandtl 
number  product  based  on  this  clearance  exceeds  about  10^,  we  should  therefore  expect 
to  begin  to  see  the  effects  of  gravity-induced  convective  currents  within  the  vessel. 

In  this  way,  the  steady  state  experiments  of  Kraussold  can  be  used  to  obtain 
a  quasi-steady  estimate  (which  is  probably  conservative)  of  the  duration  of  the  purely 
conductive  portion  of  the  cooling  curve. 

•  If  we  let  be  the  critical  value  of  t  beyond  which  bouyancy  effects  be¬ 

come  observable,  then  the  argument  above  leads  to  the  following  law 

'^crit  (19) 
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where  Gr  is  now  based  on  the  bomb  radius  R.  In  dimensional  form,  the  time  should 

then  be  independent  of  bomb  size,  being  proportional  to 

t  oc 
cnt 

or,  equivalently 

t  . .  oc  (Pn  a’/’(gpAT)"^/^  (^l) 

cnt  A 

This  last  form  can  be  used  to  estimate  an  approximate  peak  pressure  dependence, 


Thus,  the  shortest  useful  test  times  should  occur  at  the  highest  peak  pressures.  Ac¬ 
cording  to  this  model,  Lewis  and  Von  Elbe's  criterion  when  the  effect  is  most  "marked" 
correspond  to  those  cases  for  which  we  would  expect  t^^.^  to  be  large.  Vihen  the  effect 
is  not  "marked",  it  is  probably  occurring  everywhere  since  t^^.^  is  smaller  under  these 

conditions  of  conductivity  and  density. 

For  a  six  inch  diameter  bomb  giving  a  peak  pressure  of  1000  psia  in  the  oxy- 
hydrogen  system,  the  previous  estimate  leads  to  the  conclusion  that  convection  effects 
would  not  appear  before  about  10  percent  of  the  energy  had  been  lost  to  the  vessel 
walls  by  conduction,  and  this  would  take  about  20  milliseconds  to  occur.  It  should  be 
added  that  for  a  stoichiometric  hydrogen-oxygen  explosion  in  vessels  of  this  size,  the 
rise  time  (time  from  ignition  to  peak  pressure)  for  initially  atmospheric  gaseous  ex¬ 
plosions,  is  of  the  order  of  only  2  milliseconds29>3o.  For  loading  pressures  (approxi¬ 
mately  120  psia)  necessary  to  achieve  final  pressures  of  the  order  of  1000  psia,  the 
rise  time  should  be  even  less  than  2  milliseconds.  One  is  therefore  lead  to  believe 
that  diffusivity  experiments  of  the  type  investigated  here  can  be  done  without  having 
to  adopt  artificial  means  (e.g. ,  rotation)  to  suppress  or  delay  buoyancy  effects. 

Wiile  this  conclusion  is  by  no  means  certain,  it  is  probably  wisest  to  proceed  upon 
this  assumption  and  to  look  for  the  effects  of  buoyancy  in  the  longtime  pressure-records. 
Local  heat  flux  measurements  obtained  with  the  vessel  in  different  angular  positions 
could  also  be  used  to  determine  the  severity  of  these  buoyancy  effects. 


23 


8 


e.  Energy  Loss  by  Radiation  to  the  Bomb  1/fells.  In  all  gaseous  explosions, 
some  fraction  of  the  total  chemical  energy  liberated  is  radiated  to  the  walls  of  the 
containing  vessel.  In  1907,  Hopkinson*^’  observed  that  the  peak  pressure  in  coal  gas/ 
air  explosions  was  slightly  influenced  {3%)  and  the  rate  of  pressure  decay  markedly 
influenced  (505?)  when  the  walls  of  a  silver-plated  pressure  vessel  were  blackened  with 
lampback.  H.  T.  David'^^,  extending  the  work  of  Hopkinson,  continuously  monitored  the 
radiation  escaping  through  a  small  Flourite  window  in  the  bomb.  For  the  case  of  a 
25.4^  hydrogen-air  mixture  with  =  2400°K,  only  0.55?  of  the  heat  of  combustion  was 

lost  by  radiation  during  the  explosion  period  and  16.15?  during  the  entire  course  of 
the  explosion.  For  the  case  of  a  15^  coal  gas-air  mixture,  David  found  that  the  per¬ 
centage  of  the  heat  of  combustion  emitted  as  radiation  decreased  with  increased  initial 
pressure.  Recent  measurements  by  Leah  and  co-workers^^>^^  on  radiation  losses  during 
explosions  of  dry  carbon  monoxide-air  mixtures  (2  55?  CO,  755?  air)  indicate  that  at  the 
moment  of  peak  pressure  about  125?  of  the  heat  of  combustion  has  been  lost  by  radiation; 

205?  of  this  loss  is  chemiluminescent  in  nature.  For  times  of  the  order  of  (3/2)t  , 

m 

the  corresponding  values  are  17.65?  for  the  total  radiation  loss  of  which  chemilumin¬ 
escence  is  responsible  for  145?.  The  addition  of  0.52^  hydrogen  reduced  the  total  radi¬ 
ation  loss  up  to  t^  to  only  2.1^  of  the  heat  of  combustion.  Of  this,  17.85?  is  attri¬ 
buted  to  chemiluminescence.  In  the  case  of  hydrogen-air  explosions,  very  little  chemi¬ 
luminescence  has  been  observed^^.  It  appears  then  that  even  when  chemiluminescence  is 
important,  it  is  only  important  when  the  total  radiation  flux  is  small.  For  high  pres¬ 
sure  explosions,  in  which  hydrogen  is  the  combustible  gas,  the  energy  lost  to  the  walls 
by  radiation  up  to  the  moment  of  peak  pressure  probably  does  not  exceed  0.5^  of  the 
total  heat  of  combustion. 


Since  radiation  can  contribute  significantly  to  the  relative  heat  loss 

* 

Q  (t),  the  data  reduction  procedure  outlined  earlier  for  determining  the  diffusivity 

a  will  only  succeed  if  the  Rayleigh  solution  is  entered  not  with  the  total  relative 
m  ^ 

heat  loss  Q  as  determined  from  Ap(t)/Ap(oo  )  but  with  the  "conductive"  portion 
#  tt  «  * 

Q.  =  Q  "  (t)  where  (t)  represents  the  fractional  heat  loss  due  to  radiation  up 

Act  ^ 

to  time  t.  Granting  that  something  is  known  about  the  functional  form  of  ,  it  is 
possible  that  a  single  explosion  experiment  (in  which  only  the  pressure  trace  is  re¬ 
corded)  could  be  used  to  extract  not  only  the  thermal  diffusivity  of  the  gas  mixture 
at  the  moment  of  peak  pressure,  but  also  the  effective  radiative  emissivity  of  this  gas. 

Consider,  for  simplicity,  that  the  walls  of  the  enclosure  are  perfect  absor¬ 
bers.  Assume  further  that  with  the  gas  mixture  is  associated  a  total  hemispheric 
emissivity  e,  which  is  a  function  only  of  the  optical  density  of  some  species  k  and 
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the  average  temperature  in  the  vessel  (as  determined  from  the  instantaneous  pressure). 
If,  over  a  small  range  of  pressures  and  temperatures  near  p^  and  T^,  the  mole  fraction, 
X|^,  of  species  k  does  not  appreciably  change,  then  the  emissivity  e  can  be  represented 
by  a  power  law**  of  the  form 


£OC  (pR)^ 


(23) 


Since  the  radiation  heat  loss  up  to  time  t  can  be  written 

P  t  _  _ 

Q^(t)  =  4TrR2  J  e  (j  T*  dt 
0 


we  then  have. 


(24) 


»  _  4n-R^ 
Q(oo) 


m 


m 


J  (p/p„  dt 


(25) 


The  sum,  a+b,  appearing  in  the  exponent  of  p/p^  is  usually  a  small  fraction  of  4  and 
can  sometimes  be  estimated  in  advance  from  a  knowledge  of  the  principle  radiating 
species.  Wien  this  is  the  case  the  integral  appearing  in  Eg. (25)  defines  a  new 
"stretched"  time  t  v/hich  is  immediately  calculable  from  the  experimental  pressure-time 
trace.  Then  Q  will  be  linearly  proportional  to  the  product  of  t  and  the  as  yet  un- 
determined  parameter  e  .  Since  each  of  the  remaining  factors  is  calculable  from  the 

#  W  If 

known  reactants  and  observable  data,  the  quantity  Q..  computed  from  Q  -  Q  can  be 

A  E 

regarded  as  a  function  of  the  parameter  e  .  In  principle  then,  there  should  only  be 
one  value  of  e^;  i.e. ,  the  true  value  of  e^,  for  which  the  plot  of  t(Qj^  )  against 
tj^/R^  is  straight.  The  slope  of  this  line  will  be  the  true  diffusivity,  a^,  of  the 
gas  mixture  at  the  moment  of  peak  pressure.  This  is  illustrated  schematically  in- 
Fig.  14.  It  remains  to  be  seen  whether,  in  practice,  the  accuracy  of  pressure-time 
traces  will  allow  a  procedure  of  this  type  to  be  carried  to  completion. 


3.  Description  of  Apparatus 

In  view  of  the  potentialities  of  the  explosion  techniques,  it  was  decided 
to  conduct  a  series  of  preliminary  explosion  experiments  for  the  purpose  of  gaining 
familiarity  with  the  procedure  and  with  the  requisite  instrumentation.  Wiile  a  spher¬ 
ical  enclosure  would  have  been  desirable  for  direct  applicability  of  the  theory  of 
conductive  and  radiative  heat  loss,  the  more  limited  objectives  above  could  best  be 
served  by  making  use  of  a  cylindrical  vessel  which  was  available  at  the  time  and  which 


was  easily  modified  to  meet  the  experimental  requirements  of  strength  and  accessibility. 

The  apparatus  required  for  a  complete  experiment  can  best  be  broken  down, 
for  descriptive  purposes,  into  four  categories: 

1)  the  explosion  vessel  itself, 

2)  a  filling  (loading)  system, 

3)  transient  pressure  instrumentation,  and 

4)  transient  heat  flux  instrumentation. 

Several  features  of  each  of  these  systems  are  discussed  below.  This  will  be  followed 
by  a  summary  of  the  results  obtained  to  date  and  an  outline  of  possible  improvements 
in  the  experimental  technique. 

a.  The  Explosion  Vessel.  The  vessel  itself  consisted  of  a  modified  cylin¬ 
drical  reactor  of  stainless  steel  construction.  Originally,  this  vessel  had  an  inter¬ 
ior  depth  of  24  cm  and  an  internal  diameter  of  17.32  cm.  The  wall  thickness  was  3.02 
cm.  The  head  was  sealed  with  a  lead  gasket  and  contained  only  a  rupture  disc  at  its 
midpoint.  The  vescel  as  used  in  the  experiment  differed  from  the  original  in  the 
following  respects  (see  Fig.  15).  First,  the  lead  gasket  was  replaced  by  an  0-ring 
(Dash  No.  440)  seal.  The  interior  dimensions  of  the  vessel  were  altered  by  the  addi¬ 
tion  of  a  cylindrical  insert  designed  to  leave  a  cavity  whose  depth  was  equal  to  its 
diameter  (17.32  cm).  It  is  readily  verified  that  such  a  cylinder  has  the  same  surface 
to  volume  ratio  (6/d)  as  an  inscribed  sphere  of  the  same  diameter.  Fastened  to  this 
insert  was  a  0.159  cm  diameter  stainless  steel  electrode  (Ej)  which  reached  up  to  the 
center  of  the  cavity.  The  other  electrode  (E,)  consisted  of  a  similar  rod  which  ex¬ 
tended  horizontally  from  the  center  of  the  bomb  to  the  centerpost  of  an  18  mm  automobile 
spark  plug  (S)  mounted  in  the  sidewall  of  the  vessel.  Since  it  was  decided  to  use  the 
center  hole  in  the  cover  for  a  transient  heat  flux  gauge  (T)  (to  be  described),  three 
further  holes  had  to  be  machined  into  the  vessel  walls,  i .e. ,  provision  for  the  pres¬ 
sure  pickup  (P),  gas  inlet  and  exhaust  (lE)  and  condensate  drainage  (D).  The  first  two 
of  these  were  added  to  the  cover.  The  latter  was  machined  in  the  base  of  the  bomb. 
Electrically,  the  bomb  itself  was  grounded  and  the  ignition  was  accomplished  by  passing 
a  spark  across  the  centrally  located  electrode  gap.  For  this  purpose,  two  distinct 
energizing  techniques  have  been  used,  the  first  of  which  made  use  of  a  standard  auto¬ 
mobile  induction  coil.  The  second  made  use  of^a  capacitor  discharge.  In  each  case, 
completion  of  the  ignition  circuit  was  used  to  trigger  the  sweep  of  the  oscilloscope 
on  which  the  pressure  and  heat  flux  traces  were  recorded. 


b.  The  Filling  (loading)  System.  The  system  used  for  loading  and  evacuating 
the  explosion  vessel  is  shown  schematically  in  Fig.  16.  Apart  from  a  series  of  valves, 
its  main  components  are  a  pressure  gauge  (at  D),  a  vacuum  pump  (at  E)  and  the  fuel- 
oxidizer  supply  tanks  (at  C). 

The  loading  procedure  adopted  is  as  follows.  With  valves  1  through  6,  all 
open,  the  entire  system,  including  the  pressure  vessel,  is  first  evacuated  and  the  pres¬ 
sure  recorded.  The  shut-off  valve  4  is  then  closed,  isolating  the  vacuum  pump  from  the 
remainder  of  the  system.  At  this  point,  the  regulator  valve  on  gas  bottle  1  is  opened, 
and  this  gas  is  admitted  into  the  system  until  some  predetermined  pressure  level  is 
reached.  Vlhen  the  solenoid  valve  6  is  closed,  this  gas  is  then  trapped  inside  the  bomb. 
By  opening  valve  4,  the  filling  system  is  once  again  evacuated,  to  make  ready  for  the 
escape  of  gas  2.  Valve  4  is  then  closed,  and  the  regulator  valve  on  gas  bottle  2  is 
opened  to  admit  this  gas  into  the  filling  system.  The  normally  closed  solenoid  valve 
(6)  is  then  energized  and  closed  when  the  pressure  in  the  bomb  and  filling  system  is 
slightly  in  excess  of  the  desired  initial  value.  The  filling  system  is  again  evacuated 
by  reopening  valve  4.  The  metering  valve  (5)  may  then  be  used  to  bleed  off  small  am¬ 
ounts  of  gas  to  approach  the  desired  starting  pressure  p.  provided  that  enough  time  is 

t  ^ 

allowed  for  the  mixing  of  reactants  in  the  vessel  .  l/lhen  the  filling  system  is  again 
re-evacuated  or  filled  with  air,  the  bomb  is  ready  for  ignition,  subject  to  the  addi¬ 
tional  requirement.  A  knowledge  of  the  barometric  pressure,  vessel  temperature,  and 
pressures  after  the  addition  of  reactants  1  and  2,  as  well  as  after  the  bleed  period, 
enable  the  determination  of  both  the  total  mass  input  and  the  mixture  ratio.  Provision 
can  also  be  made  for  the  addition  of  diluent  gases,  if  desired. 

c.  Transient  Pressure  Instrumentation.  Commercially  available  strain  gauge 
pressure  pickups  were  used  to  obtain  the  pressure-time  history  of  the  gases  in  the 
bomb.  These  pickups  were  threaded  into  the  bomb  cover  to  a  depth  such  that  the  pres¬ 
sure  sensing  surface  was  flush  with  the  interior  of  the  vessel.  Records  have  been  ob¬ 
tained  thus  far  with  two  distinct  pickups.  The  first  of  these  was  a  Dynisco  Co. 

PT  31-2M,  0  -  2000  psi  unbonded  strain  gauge  with  a  l/4  inch  diaphragm  and  17000  cps 
resonant  frequency.  This  gauge  was  supplied  with  a  6  volt  dc  input  (excitation)  and 
had  an  experimentally  determined  sensitivity  of  98.1  psi  per  millivolt  in  the  pressure 


At  the  largest  loading  pressures  contemplated  (l25  psia),  the  calculated  mixing  time 
for  a  hydrogen-oxygen  mixture  is  about  30  minutes, based  on  molecular  diffusion  alone. 

^Dynamic  Instrument  Company,  Cambridge,  Massachusetts 


range  up  to  300  psia.  The  second  pressure  pickup  used  was  a  Dynisco  Model  PT  76B-1M 
strain  gauge  with  an  intended  range  of  0  -  1000  psi.  The  nominal  full  range  sensitivity 
of  this  gauge  was  2.99  millivolts  per  1000  psi  per  volt  excitation.  At  6  volt  excita¬ 
tion  this  corresponded  to  55.5  psi/mv.  Because  of  the  differences  in  construction  be¬ 
tween  the  aforementioned  gauges,  it  was  anticipated  that  the  second  would  be  less 
affected  by  mechanical  oscillations  in  the  vessel  walls. 

d.  Transient  Heat  Flux  Instrumentation.  If  it  is  possible  to  record  the 
temperature-time  history  T^(t)  of  a  thick  material  of  known  thermal  diffusivity,  a^, 
then  the  instantaneous  heat  transfer  rate  q"(t)to  this  surface  can  be  computed  from 
the  well-known  relation^*^ 


q”(t)  =-^\  ■  -  . . -T 

i/tt'^0  -/a  •  (t  -  t') 

Since  an  independent  knowledge  of  the  instantaneous  local  heat  flux  is  desirable,  two 
fast  response  surface  temperature  recording  methods  have  been  investigated,  with  only 
marginal  success.  The  first  of  these  {see  Fig.  17)  made  use  of  a  thin  film  (f)  of 
platinum  painted  on  a  quartz  backing  plate  (Q).  The  change  in  the  resistance  of  this 
film  was  then  used  to  deteimine  the  instantaneous  surface  temperature  of  the  quartz, 
as  outlined  in  Ref.  48  and  49.  The  second  method  made  use  of  a  commercially  available 
"surface  thermocouple"  probe^®.  This  latter  probe,  supplied  by  the  Nanmac  Corporation^, 
made  use  of  a  copper- constantan  junction  formed  across  a  thin  mica  strip  by  polishing 
the  sensing  tip.  lAhile  the  response  time  of  the  thermocouple  gauge  is  of  the  order  of 
10“^  msec,  the  uncertain  properties  of  the  backing  material  make  some  kind  of  calibra¬ 
tion  process  necessary,  if  absolute  heat  fluxes  are  desired. 

4.  Discussion  of  Preliminary  Experimental  Results 

Experience  with  the  instrumentation  described  above  has  more  than  justified 
our  initial  decision  to  postpone  the  construction  of  a  spherical  vessel  pending  the 
outcome  of  a  series  of  preliminary  experiments  with  the  existing  cylindrical  vessel. 
Perhaps  the  most  striking  feature  of  the  pressure-records,  for  stoichiometric  explosions, 
are  the  high  frequency  oscillations  ("hash")  superimposed  on  the  trace  (see  for  example, 

^Type  G  "Pencil  Probe";  Nanmac  Corporation,  P.O.  Box  8,  Indian  Head,  Maryland 


Fig.  18A).  This  particular  record  corresponds  to  the  following  initial  conditions: 

p.  =2.81  atm,  (32.9^2»  64.1^H2,  3^''N2"). 

The  gauge  sensitivity  is  55.5  psi/mv.  Oscillations  of  this  type  occurred  for  both 
pressure  pickups  over  the  pressure  range  examined  (p.  =  1  atm  -  9  atm).  The  effect 
of  diluents  was  not  examined.  Fig.  18B  shows  a  filtered  record  for  p.  =  2.95  atm,  in¬ 
itial  mixture:  yh.k%zt  65.6^2*  1^"N2".  The  gauge  sensitivity  is  98.1  psi/mv.  This 
trace  clearly  shows  the  general  shape  of  the  pressure  rise  and  subsequent  decay.  The 
horizontal  line  on  the  far  right  is  the  asymptotic  value  of  the  pressure  for  long  times. 
This  asymptote  is  lower  than  p.  due  to  the  volume  change  upon  reaction  and  condensation 
of  water  vapor  within  the  vessel.  The  quality  of  the  unfiltered  traces  obtained  thus 
far  has  not  been  adequate  for  the  data  reduction  procedure  outlined  earlier.  It  appears 
that  the  oscillations  can  be  suppressed  by  adding  diluents  to  the  initial  mixture^’. 
There  is  also  evidence  that  the  oscillation  problem  is  minimized  in  spherical  cavities. 
If  the  oscillations  occur  for  stoichiometric  mixtures  in  spherical  cavities,  then  the 
need  to  add  diluents  would  constitute  a  fundamental  limitation  of  the  experimental 
technique. 

A  large  number  of  explosions  were  carried  out  to  develop  a  thin  film  (plat¬ 
inum)  resistance  thermometer  for  making  direct  heat  flux  measurements,  Hanovia  paints 
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were  applied  to  quartz  discs  supplied  by  the  Film  Ohm  Corporation  .  These  discs  were 
then  inserted  into  a  brass  mount  (see  Fig.  17),  The  films  were  then  soldered  to  a  cen¬ 
ter-post  terminal  (T,)  and  the  mount  itself  (T2).  Records  obtained  with  these  gauges 
showed  two  types  of  defect.  First,  many  records  exhibited  a  sharp  spike  superimposed 
on  an  otherwise  reasonable  resistance-time  trace.  Secondly,  serious  zero  shifts  in 
gauge  resistance  were  observed  from  test  to  test.  Experience  to  date  suggests  that 
the  mount  shown  in  Fig.  17  is  unacceptable  in  that  the  quartz  discs  do  not  appear  to 
be  adequately  supported  and  are  thus  subjected  to  physical  strain  during  the  course  of 
the  explosion.  In  addition,  the  method  described  above  for  making  electrical  contact 
with  the  platinum  film  proved  to  be  unequal  to  the  severe  environment  to  which  the 
soldered  joints  were  exposed.  However,  if  these  mounting  difficulties  can  be  overcome, 
the  resistance  thermometer  gauge  would  have  much  to  recommend  it  over  commercially 
available  surface  thermocouple  gauges  of  the  type  described  earlier. 


Film  Ohm  Corporation,  48  W.  25th  Street,  New  York  10,  New  York 

These  discs  are  used  in  the  fabrication  of  microwave  resistors.  iVe  are  particularly 
indebted  to  Mr.  Charles  Kc  sman,  of  the  Film  Ohm  Corporation,  for  his  cooperation  and 
suggestions  with  regard  to  the  use  of  the  quartz  discs. 
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5.  Summary  of  High  Pressure  Program 


The  prospects  for  using  constant  volume  explosions  to  determine  energy  trans¬ 
port  properties  of  chemically  reacting  gases  under  conditions  of  extreme  pressures  and 
temperatures  have  been  discussed  and  possible  theoretical  and  experimental  approaches 
to  the  problem  have  been  outlined.  In  particular,  by  combining  theoretical  solutions 
for  the  energy  transfer  rate  from  a  compressible,  heat  conducting  and  radiating  sphere 
of  gas  to  the  confining  vessel  walls,  with  a  relatively  small  number  of  experimental 
measurements  (in  the  extreme  case,  transient  pressure  alone)  it  may  be  possible  to  ex¬ 
tract  both  the  thermal  diffusivity  and  radiative  emissivity  of  the  enclosed  gases.  In 
view  of  the  potentialities  of  the  technique,  a  series  of  preliminary  explosion  experi¬ 
ments  have  been  carried  out  using  stoichiometric  hydrogen-oxygen  mixtures  in  a  cylin¬ 
drical  vessel.  The  principal  difficulties  encountered  in  these  experiments  were  the 
presence  of  high  frequency  pressure  oscillations  superimposed  on  the  desired  pressure- 
time  decay  ("cooling  curve")  and  repeated  failure  of  thin  film  heat  transfer  gauges 
exposed  to  the  explosion  products.  The  pressure  oscillations  observed  in  the  present 
experiments  may  not  be  so  severe  in  spherical  cavities.  Evidence  also  exists  that  they 
can  be  suppressed  by  the  addition  of  diluents.  The  possibility  of  obtaining  reliable 
heat  flux  traces  will  then  depend  upon  improvements  in  the  mounting  of  such  gauges. 
Surface  thermocouple  gauges^^*^^,  of  the  type  used  in  automotive  practice  have  been 
investigated  and  appear  to  be  free  of  this  ruggedness  problem.  However,  such  gauges 
may  be  less  accurate  for  heat  flux  measurements  unless  they  are  carefully  calibrated 
in  a  known  heat  flux  situation. 


IV.  CONCLUDING  REMARKS 


On  the  basis  of  the  research  carried  out  thus  far,  two  potentially  fruitful 
methods  have  been  delineated  for  experimentally  studying  energy  transfer  rates  in  chem¬ 
ically  reacting*gases  under  extreme  conditions  of  pressure  and  temperature.  Prelimin¬ 
ary  results  have  been  given  both  for  chemically  reacting  ethylene/oxygen  and  hydftgen/ 
oxygen  in  the  pressure  range  5  -  15  mm  Hg  as  well  as  for  hydrogen/oxygen  explosions  at 
pressures  up  to  1000  psia.  It  is  hoped  that  the  present  work  will  stimulate  further 
investigation  along  the  lines  proposed,  with  perhaps  a  greater  awareness  of  the  poten¬ 
tial  pitfalls,  both  experimental  and  theoretical.  "...  the  subject  is  a  fair  field 
for  the  struggle  to  gain  knowledge  by  scientific  reasoning;  and,  win  or  lose,  we  find 
the  joy  of  contest"^. 
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APPENDIX  1  -  THERMAL  AND  PHYSICAL  PROPERTIES  OF  MERCURY 


Two  short  tables  of  several  useful  heat  transfer  properties  of  mercury  are 
included  here  for  convenience.  The  data  include  density,  specific  heat,  thermal  con¬ 
ductivity,  absolute  viscosity,  volumetric  heat  capacity,  thermal  diffusivity,  kinematic 
viscosity  and  Prandtl  number  for  heat  conduction  over  the  temperature  range  0  to  250°C, 
at  atmospheric,  or  near-atmospheric,  pressure. 

TABLE  IV 

HEAT  TRANSFER  PROPERTIES*  OF  MERCURY 

Density,  Specific  Heat,  Thermal  Conductivity 
and  Absolute  Viscosity  at  Total  Pressures  Near  Atmospheric 


T°C 

P 

c 

P 

100  X 

t\ 

0 

13.60 

0.03346 

1.96 

1.68 

273 

20 

13.55 

0.03325 

2.09 

1.55 

293 

40 

13.50 

0.03308 

2.20 

1.44 

313 

60 

13.45 

0.03294 

2.31 

1.35 

333 

80 

13.40 

0.03284 

2.42 

1.27 

353 

100 

13.35 

0.03269 

2.51 

1.21 

373 

200 

13.11 

0.03230 

2.96 

1.01 

473 

250 

13.00 

0.03210 

•  •  •  » 

0.96 

523 

^Liquid  Metals  Handbook,  Second  Edition,  Government  Printing  Office,  Washington,  D.C., 
1952;  see  also  Handbook  of  Chemistry  and  Physics,  Chemical  Rubber  Publishing  Company, 
Cleveland,  Ohio 


Units  are  as 
conductivity 


follows:  density  in  gm(cm)"^;  specific  heat  in  cal(gm)"’ >  thermal 
in  (cal) (cm)"^(sec)"M°K/cm)"’ ;  absolute  viscosity  in  gm(cm)"Msec)"’ 


@ 


36 


* 


9 


TABLE  V 


HEAT  TRANSFER  PROPERTIES^  OF  MERCURY 

Volumetric  Heat  Capacity,  Thermal  Diffusivity,  Kinematic  Viscosity 
and  Prandtl  Number  at  Total  Pressures  Near  Atmospheric 


T°C 

pc 
^  P 

a 

V 

(Pr,F 

T°K 

0 

0.4549 

0.0431 

0.00124 

0.0237 

0.169 

273 

20 

0.4504 

0.0464 

0.00114 

0.0246 

0.157 

293 

40 

0.4465 

0.0493 

0.00107 

0.0216 

0.147 

313 

60 

0.4430 

0.0522 

0.00100 

0.0192 

0.139 

333 

80 

0.4401 

0.0549 

0.000949 

0.0173 

0.132 

353 

100 

0.4365 

0.0576 

0.000906 

0.0157 

0.125 

373 

200 

0.4236 

0.0698 

0.000770 

0.0110 

0,105 

473 

250 

0.4172 

•  •  • 

0.000739 

»  •  • 

•  •  • 

523 

Units  are  as  follows:  volumetric  heat  capacity  in  cal(cm)“^(°K)"’ ;  thermal  diffusivity 
and  kinematic  viscosity  in  {cm)^(sec)"’ ;  Prandtl  number,  dimensionless 


APPENDIX  2  -  EFFECTS  OF  WATER  VAPOR  CONDENSATION  ON  THE  EXPLOSION  METHOD 


Condensation  of  water  vapor  is  one  of  the  causes  of  variations  in  the  spec- 
ific  heat,  c^,  during  the  cooling  process.  Moreover,  the  diffusion  of  water  vapor  to 
the  cold  wall,  where  it  condenses  exothermically,  constitutes  a  potentially  important 
energy  transfer  mechanism  other  than  ordinary  heat  conduction.  This  results  in  the 
well-known  calorimetric  concept  of  "higher"  and  "lower"  heating  values  for  fuels  con¬ 
taining  hydrogen,*  i  .e. ,  each  fuel  will  have  two  heating  values,  which  will  differ  by 
the  product  of  the  latent  heat  of  vaporization,  H,  and  the  total  mass  of  water  formed 
by  condensation  in  the  vessel.  With  regard  to  use  of  an  explosion  method  for  deducing 
gas  phase  thermal  diffusivities,  it  is  intuitively  clear  that  if  the  difference  between 
the  higher  and  lower  heating  values  of  tte  combustible  is  negligible  compared  to  the 
absolute  value  of  the  lower  heating  value,  then  the  effect  of  water  vapor  condensation 
on  the  inferred  thermal  diffusivity  will  necessarily  be  small.  If  the  condensation 
occurs  mainly  at  the  walls  themselves,  then  the  fact  that  the  Fick  diffusivity 

is  generally  less  than  the  thermal  diffusivity  a  of  the  gas  mixture  would 

n^U^rnlX 

further  contribute  to  the  relative  importance  of  energy  transport  by  ordinary  heat 
conduction.  This  is  particularly  clear  in  the  case  of  diffusion  controlled  vaporiza¬ 
tion  for  short  times,  since  the  square  root  law  immediately  leads  to  the  con¬ 

clusion  that,  up  to  any  time  t,  the  ratio  of  the  total  heat  Qjj  (liberated  as  a  result 
of  water  vapor  diffusion  and  condensation)  to  Qj^  is  given  by 

Vx  ■* 

For  longer  times,  droplets  probably  form  on  the  wall.  These  droplets  would  increase 
in  size  until  gravity  would  force  them  to  run  down  the  wall.  On  clean  surfaces  con¬ 
tinuous  films  are  also  possible.  In  both  cases,  the  presence  of  water  at  the  surface 
introduces  an  additional  resistance  for  the  flow  of  heat  to  the  wall.  Only  in  certain 
extremely  simple  cases  are  the  laws  governing  heat  transfer  in  this  type  of  environment 
theoretically  understood  (see,  for  example.  Chapter  12  of  Eckert  and  Drake"^®).  It  is 
also  clear  that  the  effective  radiative  absorptivity  of  the  vessel  walls  is  likely  to 
change  as  condensation  proceeds.  For  these  reasons,  unless  this  process  can  bo  arti¬ 
ficially  inhibited  (delayed),  the  accuracy  of  the  explosion  method  as  a  means  of  deter¬ 
mining  gas  phase  transport  properties  is  expected  to  degenerate  in  the  presence  of  large 
amounts  of  condensation. 
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FIG.  1  CONFIGURATION  OF  BURNER  AND  HEAT  TRANSFER  PROBE  (SCHEMATIC) 


FIG.  3  CFSDSS-SECTIONAL  DETAIL  OF  PROBE  THERMOCOUPLE  AND  COOLANT  FLOW 
CONFIGURATION  (TO  SCALE) 


CONSTANT  HEAD  MERCURY  FLOW  SYSTEM  (SCHEMATIC) 


FIG.  5  SCHEMATIC  OF  LOW  PRESSURE  BURNER  ASSEMBLY 


LUMINOUS  ZONE 


DISTANCE  FROM  SCREEN  TO  PROBE  CENTERLINE,  cm 
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FIG.  9  COOLING  PROCESS  REPRESENTED  ON  THE  INTETOIAL  ENERGY-DENSITY  PLANE 


C  OXY-HYDHOGEM  EXPLOSIONS 
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FIG.  11  DEPENDEFCE  OF  MOBMALIZED  HEAT  LOSS  ON  NDNDIMENSIONAL  TIME 

(FOURIER  QUENCHED  SPHERE  SOLUTION) 


graphical  determination  of  effective  thermal  diffusivity 

FFDM  RELATIVE  CONDUCTIVE  HEAT  LOSS  MEASUREMENTS 


FREE  CONVECTION  IN  ENCLOSED  SPACES 

STEADY-STATE,  WITH  INNER  BOUNDARY 
TRANSIENT,  WITH  NO  INNER  BOUNDARY 


FIG.  14  GRAPHICAL  DETERMINATION  OF  BOTH  DIFFUSIVITY  AND  EMISSIVITY 

FROM  PRESSURE  MEASUREMENTS  ALONE 


FIG.  16  EVACUATION  AND  LOADING  SYSTEM 


FIG.  17  CRDSS-SECTION  OF  THIN  FILM  RESISTANCE  THERMOMETER  ASSEMBLY 
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TIME*  msec 


A.  Initial  Pressure  =  2.81  atmi  Initial  Composition  =  32.95f  0*,  64.1$  Hj,  3$  "Nj"| 
Gauge  Sensitivity  =  55.5  psi/mv 


B.  Initial  Pressure  =  2.95  atmi  Initial  Composition  =  33.4$  0„  65.6$  H,,  1$  '*Na"| 
Gauge  Sensitivity  =98.1  psi/mv 


FIG.  18  PRESSURE  PICKUP  OUTPUTS  AS  A  FUNCTION  OF  TIME 
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